Temperature dependence of line parameters of 12C16O2 near 2.004 μm studied by tunable diode laser spectroscopy
Ma Hongliang1, 2, Sun Mingguo1, Zha Shenlong2, Liu Qiang1, Cao Zhensong1, †, Huang Yinbo1, Zhu Zhu2, Rao Ruizhong1
Key Laboratory of Atmospheric Optics, Chinese Academy Sciences, Hefei 230031, China
Anqing Normal University, Anqing 246133, China

 

† Corresponding author. E-mail: zscao@aiofm.ac.cn

Abstract

The absorption spectrum of carbon dioxide at has been recorded at sample temperatures between 218.0 K and room temperature, by using a high-resolution tunable diode laser absorption spectrometer (TDLAS) combined with a temperature controlled cryogenically cooled absorption cell. The self-, N2-, and air-broadening coefficients for nine transitions of belonging to the band in the 4987 cm−1–4998 cm−1 region have been measured at different temperatures. From these measurements, we have further determined the temperature dependence exponents of the pressure-broadening coefficients. To the best of our knowledge, the temperature dependence parameters of the collisional broadening coefficients are reported experimentally for the first time for these nine transitions. The measured halfwidth coefficients and the air temperature dependence exponents of these transitions are compared with the available values reported in the literature and HITRAN 2012 database. Agreements and discrepancies are also discussed.

1. Introduction

Carbon dioxide (CO2) is the major anthropogenic greenhouse gas that dramatically increased since the beginning of the industrialization era due to the effects of human activities on the Earth’s environment.[1] Its infrared absorption and radiation play a crucial role in radiative transfer processes especially for the evolution of the Earth’s atmosphere and climate. Therefore, accurate spectroscopic data, in particular in the near-infrared region, is essential, e.g., for reliable quantitative studies of global warming or analysis of remote sensing data.[2]

Currently, the global monitoring of CO2 is of high importance and several methods[37] are developed for the monitoring of atmospheric carbon dioxide. Particularly, in 2009, the Greenhouse Gases Observing Satellite (GOSAT)[8] was launched by Japan to provide global time-dependent maps of CO2 concentration in the 1.6- and 2.0- region. Therefore, the determination of the exact values of the spectroscopic parameters of the spectral lines around is very important for the high precision measurements of CO2 concentration in atmosphere, especially the temperature dependence exponents because of the wide temperature range in the Earth’s atmosphere. Due to the importance in remote sensing applications, the absorption spectrum of CO2 around at different temperatures has been quantitatively studied by many groups.[914] In 1979, Valero et al.[9] reported the line intensities and half-widths of 51 lines at different temperatures for the band of CO2 at 4854 cm−1 with a Fourier transform spectrometer (FTS) with resolution of 0.075 cm−1. Subsequently, in 1990, Suárez et al.[10] published the temperature dependence of self-broadened half-widths of CO2 with an FTS and at a resolution of 0.004 cm−1 in the range 4850 cm−1–5050 cm−1. Recently, Joly L et al.[11] and Li J S et al.[12] determined the line strength, air-broadening half-width, air pressure-induced shift coefficient as well as their temperature dependences in the range 4870 cm−1–4880 cm−1 of CO2, in order to provide reference data for the differential absorption lidar instrument. Li J S et al.[13] also reported the self-induced pressure shift coefficients and the temperature dependence for 8 transitions in this region by a high resolution tunable diode laser absorption spectrometer combined with a cryogenically cooled optical multi-pass cell in 2012. More recently, the CO2 line parameters including temperature dependences of air- and self-broadened line shapes near 2.06 have been studied by Benner et al.[14] Besides these results, no further low temperature experimental information has been obtained for these lines. However, the intercomparison data in Ref. [15] (Table 1, Table 3, and Table 4) show more measurements are needed for the CO2–N2 and CO2–CO2 systems, especially the temperature dependence exponents of CO2 transitions near belonging to the band, and hence the present measurements are important.

In the present work, we have investigated nine lines of CO2 in the 4987 cm−1–4998 cm−1 spectral region belonging to the band, using a narrow linewidth tunable diode laser combined with a temperature controlled cryogenically cooled absorption cell. The temperature dependence exponent n and the pressure-broadened half-width coefficients for nine transitions of CO2 with intensities stronger than ) have been determined and compared with the available values reported in the literature and the HITRAN2012[16] (hereafter shortened to HITRAN) database.

2. Experiment details

The CO2 molecular spectra in the 2.004- region are recorded with a high resolution tunable diode laser absorption spectrometer. The experimental setup is schematically shown in Fig. 1. A commercial distributed feedback (DFB) diode laser purchased from Nanoplus has been used as the light source. The side mode suppression ratio is higher than 40 dB. The diode laser emits radiation near with no mode-hops over the tunable range. The typical output power is about 7.8 mW, with a relatively narrow linewidth of about 3 MHz, which can be neglected comparing to Doppler limited molecular absorption linewidth (for the CO2 lines around at 296 K, the typical Doppler half-width at half-maximum is about 230 MHz). The continuous tuning range (at constant temperature) is within 5 cm−1. This point is of particular interest to reconstruct properly the baseline, as the laser sweeps the molecular transition over a spectral range that is large enough to yield zero-absorption signals at the end and at the beginning of the scanning interval. The laser wavelength is temperature-stabilized by means of a Peltier thermo-element and is driven by a low noise current supply (ILX Lightwave LDC-3724C). The controller was linked to a computer via a GPIB interface while its current was scanned in a step-by-step mode with typical size of 0.1 mA in order to scan the laser over the selected absorption lines by modulation of the driving current. The laser beam is usually collected by a fixed focus collimator and is separated into three parts via a wedged CaF2 beam splitter. The main beam was passed through a sample cell for the absorption signal. The second beam was introduced onto a wavemeter (Bristol-621A) to record frequency of the absorption signal with the accuracy of 0.001 cm−1 at , and the third one was used as a background signal. Two InGaAs amplified detectors (Thorlabs, PDA10DT-EC) were used to record the sample and the background spectra. A three-channel acquisition system ensures simultaneous recording of the gas sample absorption spectrum and the other two signals (wavemeter and background signal).

Fig. 1. (color online) The schematic diagram of the experimental setup for the CO2 spectra measurement.

A home-made cooling cell was developed for the purpose of atmospheric applications, and more precisely described in recent publications.[17,18] Briefly, the absorption path length is 17.4 cm and the cooling system is realized by a dry nitrogen flow, circulating around the absorption cell. Two PT100 sensors mounted at different locations on the wall of the sample cell were used to measure the gas temperature during the recordings. The gas pressure is continuously measured by two capacitance gauges (Shanghai Zhentai CPCA-200, full range 20 kPa; MKS 626C, full range 67 kPa), and the pressures are measured with an uncertainty of 0.25%. The gas samples (carbon dioxide, air, nitrogen) were bought from Nanjing Special Gas Co., with a nominal purity of 99.99%.

During the measurements, the sampling cell was first evacuated to the order of 10−2 Pa using a mechanical together with a molecular pump (Leybold, TW300), then filled with the natural abundance CO2 gas of 99.99% purity or the mixed gas (dry air+CO2 or N2+CO2). In order to mix the gases fully, absorption signals were always recorded after at least 3 hours of preparing the gas mixtures. All measurements were achieved with a temperature varying less than 0.3 K during a spectrum acquisition. For a given couple of pressure and temperature, 3–5 spectra were recorded to reduce systematic errors. Table 1 summarizes the experimental conditions for the entire set of spectra used in this work.

Table 1.

Summary of experimental conditions of CO2 spectra analyzed in this work.

.
3. Data processing

For a monochromatic laser beam, the transmission of light radiation Tν at a frequency ν through a uniform gas can be given by the Lambert–Beer law:

where I0 and are referred to the incident and transmitted power, respectively, (in unit cm−1) is the spectral absorption coefficient; L (in unit cm) is the optical path length within the cell.

To retrieve the Lorentz half-width coefficients of these transitions of CO2, the local baseline of the observed spectra was assumed to be a third order polynomial of the wavenumber and spectral lines were fitted with a Voigt function having four parameters, (line position, integrated absorbance, FWHM of the Gaussian and Lorentz components) that were determined using a multi-peak fitting program (http://www.unipress.waw.pl/fityk/). This software is an interactive nonlinear least-squares fitting program based on the Levenberg–Marquardt algorithm that could be simultaneously fitting a few lines for one spectrum at a time. During the fitting procedure, the FWHM of the Gaussian component of the lines were fixed to their theoretical values. The integrated absorbance were also fixed to the values that were calculated from line intensities (note that these line intensities were taken from HITRAN). Note that other weaker lines, belonging to the band of CO2, which need to be included in the nonlinear least-squares multi-peak fitting procedure for a correct determination of the line parameters. The four Voigt function parameters of these weaker lines were calculated from HITRAN. Figure 2 shows an example of the recorded spectrum and its corresponding Voigt fit.

Fig. 2. (color online) Example of recorded absorption spectrum for the R(20) transition of CO2 perturbed by air. Experimental profile and fitted Voigt profile are presented here. The observed minus calculated residuals are shown at the bottom. Pressure is 30.944 kPa and temperature is 253.2±0.2 K.

The self-broadened half-widths were fitted assuming a linear pressure dependence using the relation:

In Eq. (2) is the retrieved broadening coefficient at 1 atm ( ), is the measured broadened half-width of the spectral line and P is the total sample pressure.

For a given perturber, the collision half-width can be defined by:

where, is the partial pressure of CO2, is the coefficient for self-broadening, is the perturber (air or nitrogen) pressure, is the coefficient for collisional broadening by air and nitrogen.

Using Eqs. (2)–(3), the broadening coefficients , , and are determined by plotting the Lorentz widths versus gas pressure and then fitting them by straight lines. The collisional broadening coefficient at a given temperature is expressed as the slope of the best fit line. Figure 3 illustrates the plot of the measured air-broadened half-widths versus pressure at 218.4 K for the R(20) line in the band of CO2. The same procedure is also followed for the other transitions.

Fig. 3. (color online) Plot of collisional half-width at half maximum versus each pressure for the R(20) transition of CO2. The slope of the regression line corresponds to the air-broadening coefficient in at 218.4 K.

The temperature dependence of the half-width is expressed using the power law model, and is generally expressed as:

where is the broadening coefficient of the line at 1-atm pressure and reference temperature T = 296 K, and is the retrieved broadening coefficient at 1atm pressure and temperature T; n is the temperature dependence exponent. Figure 4 shows the air-broadening coefficient temperature dependence exponent of the R(28) transition of CO2 as an example.

Fig. 4. (color online) Fits of as a function of for the R(28) line of CO2 using air-broadened widths measured over the range 218.4 K–295.9 K.
4. Results and discussion

Within the tunable range of the Nanoplus diode laser, nine strong lines of CO2 in the band lying between 4987 cm−1 and 4998 cm−1 were investigated in the laboratory.

4.1. Self-, N2-, and air-broadening coefficients

In Fig. 5 the broadening parameters obtained with the Voigt profile are plotted as a function of the rotational quantum number m (m = J + 1 for R branch lines) at various temperatures for the R(12) to R(28) of band of CO2. In Fig. 5, we have noticed that: firstly, all the pressure-broadening coefficients increase with the diminution of the temperature; secondly, at the same temperature was noticed from the dataset; third, all the results present similar m-dependences. The error bars are one standard deviation given by the linear regression. However, these values are not considered to represent the absolute uncertainties because there may be systematic uncertainties arising from additional sources, such as those related to the knowledge of wavenumber calibration errors, gas sample pressures and normalization uncertainties. The absolute accuracies in our measured parameters due to various errors are estimated to be better than 4%, as stated from our experience.

Fig. 5. (color online) Self-, N2-, and air-broadening coefficients as a function of the rotational quantum number m at various temperatures for R(12) to R(28) of band of CO2.

In order to assess the quality of our retrieved results, we compared our broadening coefficients at room temperature with other corresponding measurements[16,1922] for this band of CO2 transitions, shown in Fig. 6(a). To extend our comparison in Fig. 6(b) we have plotted against m the percent differences {(Present work–other results)/other results×100%} between our values and other reported results. Firstly we compare our self-broadening coefficients with previous determinations of Régalia-Jarlot et al.,[19] Corsia et al.,[20] and HITRAN.[16] As shown in Fig. 6, our values are close to those of Régalia-Jarlot et al.,[19] the difference lies within −0.76% to 1.19% and a value of is achieved for the RMSD (root-mean-square error) of our results and the measurements reported by Régalia-Jarlot et al.[19] We can see that the difference between our self-broadening coefficients and the values from the HITRAN database varies from −3.7% to −6.0% for all the lines (for these transitions, the HITRAN database has assumed the self-broadening coefficients from Ref. [21]). However, the differences between our results and HITRAN (or Ref. [21]) maybe attributed to the different setup and the analysis techniques employed in the two studies. It can be noticed that in Ref. [21], the self-broadening coefficients were obtained using an FTS and the spectral resolutions for their spectra are 0.010 cm−1 to 0.013 cm−1 while we were using a very narrow laser source with a resolution of about 0.0001 cm−1. When comparing our results (for R(22), R(24), R(26), R(28)) with those from Corsia et al.[20] we found that the difference becomes large (up to −8.66%), namely, our results are systematically lower than those obtained by Corsia et al.[20] However, we note that the data from Corsia et al.[20] are clearly higher for the transitions (R(22), R(24), R(26), R(28)) in comparison with all results. At present, we cannot explain these different findings, maybe due to unknown reasons.

Fig. 6. (color online) (a) Comparison between our values of pressure-broadening coefficients in the band of CO2 obtained at room temperature and previous results obtained by other authors[19,20,22] and HITRAN.[16] (b) The percent differences {(Present work–other results)/other results ×100%} of pressure-broadening coefficients between this work and the results from Refs. [19], [20], [22], and HITRAN[16] for transitions in band are plotted as a function of the m value (Note: PW = present work).

As to N2-broadening coefficients, we could find only the values listed by Corsia et al.[20] for direct comparison. Only four lines are identical to the two works. For the transitions R(24), R(26), and R(28), the differences are −3.58%, −0.55%, and −3.14%, respectively. Regarding the R(22) line, the discrepancy in N2-broadening coefficient between our measurement and the result from Corsia et al.[20] is up to 9.36%. It is difficult to explain this rather large difference for the transition. However, it is interesting to note that, the self-broadening coefficient and O2-broadening coefficient of the R(22) line are always larger than the other three lines (R(24), R(26), and R(28)) as shown in Ref. [20].

For air-broadening coefficients, we have also compared our values with those from HITRAN and Toth et al.[22] Our results of measured air-broadening coefficients are much closer to those of HITRAN and the results of Toth et al.[22] The agreement between our values and those obtained by Toth et al.[22] is better than 3.4% for all transitions. (Note: the air-broadening coefficient of R(14) was not given in Ref. [22].) With the exception of R(12), our broadening parameters agree within 4.0% with those from HITRAN. The RMSD of our results and the measurements reported by Toth et al.[22] and HITRAN are and , respectively.

4.2. Temperature dependence exponent

For remote sensing applications, both temperature and pressure dependences of half-width coefficients are required on a line-by-line basis. In determining the temperature dependence exponents for the Lorentz self-, air-, and N2-broadened half-width coefficients, we have employed the power law relationship given in Eq. (4) by fitting the pressure-broadening coefficients at all sample temperatures. The results are summarized in Table 2. The uncertainties in parentheses represent only the uncertainties as obtained from the linear regression. The absolute uncertainties of the temperature dependence exponents are difficult to assess, but are less than 10%, as considering all the unknown sources of errors. From Table 2, unlike broadening coefficients, we can see that the temperature dependence exponents are less dependent on m values. The mean values of the temperature dependence exponent for self-, air-, and N2-broadening coefficients are 0.685, 0.660, and 0.714, respectively. Meanwhile, we also compared our air temperature dependence exponents with the values from HITRAN, as it is shown in Table 2. The RMSD of our air temperature dependence exponents and the results given in HITRAN is 0.0153, while the agreement is better than 2% for most of the transitions. Although the range in and is generally close, there is less scatter in the results of HITRAN. Yet, the dependence in the rotational quantum number m of HITRAN is a theoretical calculation, while our experimental results were affected by other factors. Unfortunately, we did not find the self (or N2) temperature dependence exponents for comparison.

Table 2.

Self-, air-, and N2-temperature dependence exponents for the R(12) to R(28) of band of CO2. The air temperature dependence exponent from HITRAN also presented in the table and compared with those of our results. Numbers in parentheses of our work are the error obtained on the linear regression.

.
5. Summary

We measured nine CO2 self-, N2-, and air-broadening half-width coefficients for the band in the 4987 cm−1–4998 cm−1 region using a narrow linewidth diode laser associated with a home-made cryogenic cell. The results show that the broadening coefficients present clear temperature- and m-dependences. Furthermore, the parameters obtained at room temperature were thoroughly compared to former studies and the HITRAN database, and some agreements and discrepancies have been found. Self-broadening coefficients were determined for nine transitions and found to be from 3.7% to 6.0% lower than those from HITRAN, but very close to those obtained by Régalia-Jarlot et al.[19] For N2-broadening coefficients, only 4 lines are available in Ref. [20] for a direct comparison. The results show that good agreement was achieve between our values and those from Corsia et al.[20] except R(22). We have also compared our air-broadening coefficients with those from HITRAN and Toth et al.[22] The agreement is within 3.4% (nine lines for Toth et al.[22]) and 3.5% (except R(12) and R(14) lines for HITRAN) for most of the transitions, respectively. From these broadening coefficients at different temperatures, the temperature dependence exponents have been determined and no significant tendency was found. Our air temperature dependence exponents agree very well with those in HITRAN with an RMSD of only 0.0153. To the best of our knowledge, the temperature dependence exponents of the self-, N2-, and air-broadening coefficients are reported experimentally for the first time for these nine transitions of CO2. We believe that the reported data set presented here will be helpful for the Earth’s radiation balance or remote sensing of CO2 in the atmospheric spectral window at .

Reference
[1] IPCC. Climate Change 2007: the physical science basis. Contribution of working group I to the fourth assessment report of the intergovernmental panel on climate change; 2007.
[2] Rothman L S Jacquemart D Barbe A 2005 J. Quantum Spectrosc. Radiat. Transfer 96 139
[3] Durry G Amarouche N Zéninari V Parvitte B Lebarbu T Ovarlez J 2004 Spectrochim. Acta 60 3371
[4] Pouchet I Zéninari V Parvitte B Durry G 2004 J. Quantum Spectrosc. Radiat. Transfer 83 619
[5] Dufour E Bréon F M Peylin P 2004 J. Geophys. Res. 109 D09304
[6] Tolton B T Plouffe D 2001 Appl. Opt. 40 1305
[7] Kuang Z Margolis J Toon G Crisp D Yung Y 2002 Geophys. Res. Lett. 29 1716
[8] Inoue G Yokota T Oguma H Higurashi A Morino I Aoki T 2004 AGU 2004 fall meeting December 13–17, 2004 San Francisco, California, USA
[9] Valero F P J Suárez C B Boese R W 1979 J. Quantum Spectrosc. Radiat. Transfer 22 93
[10] Suárez C B Valero F P J 1990 J. Quantum Spectrosc. Radiat. Transfer 43 327
[11] Joly L Marnas F Gibert F Bruneau D Grouiez B Flamant P H Durry G Dumelie N Parvitte B Zéninari V 2009 Appl. Opt. 48 5475
[12] Li J S Durry G Cousin J Joly L Parvitte B Flamant P H Gibert F Zéninari V 2011 J. Quantum Spectrosc. Radiat. Transfer 112 1411
[13] Li J S Durry G Cousin J Joly L Parvitte B Zéninari V 2012 Spectrochimica Acta Part. 85 74
[14] Benner D C Devi V M Sung K Brown L R Miller C E Payne V H Drouin B J Yu S S Crawford T J Mantz A W Smith M A H Gamache R R 2016 J. Mol. Spectrosc. 326 21
[15] Gamache R R Lamouroux J Blot-Lafon V Lopes E 2014 J. Quantum Spectrosc. Radiat. Transfer 135 30
[16] Rothman L S Gordon I E Babikov Y L 2013 J. Quantum Spectrosc. Radiat. Transfer 130 4
[17] Ma H L Liu Q Cao Z S Chen W D Vicet A Huang Y B Zhu W Y Gao X M Rao R Z 2016 J. Quantum Spectrosc. Radiat. Transfer 171 50
[18] Ma H L Sun M G Liu A W Vicet A Chen W D Cao Z S Wang G S Liu Q Gao X M Rao R Z 2015 Acta Phys. Sin. 64 153301 (in Chinese)
[19] Régalia-Jarlot L Zéninari V Parvittea B Grossela A Thomas X Heydena P von der Durry G 2006 J. Quantum Spectrosc. Radiat. Transfer 101 325
[20] Corsia C D’Amato F DeRosa M D Modugno G 1999 Eur. Phys. J. 6 327
[21] Toth R A Brown L R Miller C E Devi V M Benner D C 2006 J. Mol. Spectrosc. 239 243
[22] Toth R A Brown L R Miller C E Devi V M Benner D C 2007 J. Mol. Spectrosc. 242 131